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We apply the perturbation theory integral representation (PTIR) to solve for 
the bound state Bethe-Salpeter (BS) vertex for an arbitrary scattering kernel, 
without the need for any Wick rotation. The results derived are applicable to 



Q I any scalar field theory (without derivative coupling). It is shown that solving 

' directly for the BS vertex, rather than the BS amplitude, has several major 

advantages, notably its relative simplicity and superior numerical accuracy. In 
^ . order to illustrate the generality of the approach we obtain numerical solutions 

using this formalism for a number of scattering kernels, including cases where 
the Wick rotation is not possible. 
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I. INTRODUCTION 



We present here an improved approach to the solution of the scalar-scalar Bethe-Salpeter 
(BS) equation directly in Minkowski space, utilising the Perturbation Theory Integral Rep- 
resentation (PTIR) of Nakanishi . The PTIR is a generalised spectral representation for 
n-point Green's functions in Quantum Field Theory. 

This work extends and improves earlier work which applied the PTIR approach to the 
BS amplitude 0. Here we formulate a real integral equation for the BS vertex. This 
considerably simplifies the expression for the kernel function relative to those obtained for 
the BS amplitude [Q. In particular, some singularity structures which were present in the 
kernel of the BS amplitude equation due to the constituent particle propagators are absent 
in the corresponding expression for the BS vertex. Consequently, it is much simpler to 
implement the problem numerically, and we therefore do not encounter previous difficulties 
with residual numerical noise. We have checked that our numerical results are in agreement 
with those obtained in Euclidean space by other authors (for example. Linden and Mitter 0, 
or more recently, Nieuwenhuis and Tjon ||TO| , p!T| ). In particular, with sufficient computer time 
we have seen no limit to the accuracy that can be achieved with our formalism. We can 
routinely obtain four figure accuracy on a workstation. 

In this work we will deal exclusively with scalar theories. For simplicity we will consider 
here bound states with equal-mass constituents, although it is easy, and desirable in many 
applications, to generalise the approach to unequal-mass constituents. 

We illustrate the BS equation for a scalar theory in Fig. |I], where P) is the BS 
amplitude, P = pi + p2 is the total four-momentum of the bound state and p = rj2Pi — ?7iP2 
is the relative four-momentum for the two scalar constituents. We have then M = a/P^ fQj- 
the bound state mass and also rji + t]2 = 1, but otherwise the choice of the two positive real 
numbers t]i and 772 is arbitrary. As in Ref. [|l|] we choose here rji = 772=1/2. 

The renormalised constituent scalar propagators are D{pl 2) and K{p, q; P) is the renor- 
malised scattering kernel. For example, in simple ladder approximation in a (p'^a model we 
would have K{p, q; P) = {ig){iD„{\p — qY)){ig), where D^ijP') = — + ie) and m„ is 
the (T-particle mass. Note that the corresponding proper (i.e, one-particle irreducible) vertex 
for the bound state is related to the BS amplitude by $ = {iD){iV){iD). 

We follow standard conventions in our definitions of quantities, (see Ref. and also, 
e.g., 0). Thus the BS equation for any scalar theory can be written as 

^r(pi,P2) = / m<luq2)) mql)) K{p,q;P), (1) 

where similarly to pi and p2 we have defined qi = rjiP + q and q2 = 772 -P — 1- Equivalently 
in terms of the BS amplitude we can write 

C(pf)-' *(p, P) C(pi)-' = - / ^ *(9. P) K(P, 9; P) (2) 

where the kernel function defined by /(p, g; P) = —iK{p, q; P) is the form typically used by 
Nakanishi [Q. In ladder approximation for a (j)^a model we see for example that /(p, g; P) = 



2 



g'^ / (m^ ~ {p — (lY — T'^) ■ In this treatment we will solve the vertex version of the BS equation, 
i.e. Eq. ([l|), for an arbitrary scattering kernel. 



II. PTIR FOR SCALAR THEORIES 

In contrast to the approach in Ref. we will begin with the Bethe-Salpeter equation 
for the bound-state vertex, Eq. (|I|), rather than the equivalent equation for the amplitude, 
Eq. (^, and derive a real integral equation for the BS vertex. To do this we require spectral 
representations for the vertex, for the propagator, and for the scattering kernel / of Eq. 
(0). The renormalised propagator may be written as 

Z,(,) = -(^J^+r ,la^*^], (4) 

\m^ — q^ — ie J{m+fj)'^ a ~ — J 

where P(f,{a) is the renormalised spectral function. Note that p^{a) > 0, (see, e.g., Ref |^). 

The Bethe-Salpeter amplitude <I>(p, P) for the bound state of two ^-particles having the 
total momentum P = pi + p2 and relative momentum p = {r]2Pi — i]iP2) can be defined as 

{0mx^)cPix2)\P) = e-^^-^(O|T0(r/2a:)0(-r/ix)|P) = e^^'^ | _^e-^-$(p,P) , (5) 

where the fields for the scalar constituents are denoted by (p, and where we have made use 
of the translational invariance of the BS amplitude. Following the conventions of Itzykson 
and Zuber |Q (e.g., pp. 481-487), we define centre-of-momentum and relative coordinates 
X = TjiXi + 7/2X2 and x = Xi — X2 such that Xi = X + ri2X, X2 = X — rjix, and P ■ X +p ■ x = 

Pl-Xi+P2- X2. 

Equivalently to Eq.(|), we can write 

$(p,P) = e'^-^y d''xe''P-%Q\T(t){xi)(t){x2)\P) = J d^xe'P-^0\T(l){r]2x)(j){-riix)\P) . (6) 

Note that the bound states are normalised such that (P|P') = 2ujp{2ttY6^{P' — P), where 
ujp = (P^ + M^)^/^ with M the bound state mass. For a positive energy bound state we 
must have P^ = M^, < P^ < (2m)^, and P° > 0. The normalization condition for the BS 
amplitude is given by 

Iwi-J ^,^q.P)^{D-\pl)D-\pl){2'aY6\p-q) + K{p,q-P)] 

(7) 

where the conjugate BS amplitude <l>(p, P) is defined by 

$(p,P) = e-'^-^y d''xe-'P-%P\T(t)\xi)(t)\x2)\Q) = J d^xe-''' ^P\T(P\ri2x)(l)\-riix)\0) . 

(8) 
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A. PTIR for Scattering Kernel 



The scattering kernel I{p, q; P) = —iK{p, q; P) describes the process 00 00, where p 
and q are the initial and final relative momenta respectively. It is given by the infinite series 
of Feynman diagrams which are two-particle irreducible with respect to the initial and final 
pairs of constituent particles. For purely scalar theories without derivative coupling we 
have the formal expression for the full renormalised scattering kernel M 



oo 



n 



1 - 



— le 



1 - 



— le 



(9) 



where q^ is the 4-momentum squared carried by 0j and s,t and u are the usual Mandelstam 
variables. The symbol VL denotes the integral region of such that = {^i | < < 
1) 'Ci = 1 = 1; • • • ) 6)}. The scattering kernel PTIR can be re-written in a more compact 
form as 

/(p, q;P)=T.JdlJ di _ rS}^d% + e.,a.P + Uv.P) - le ' 



ch 



where the subscript "ch" indicates which channel we are dealing with (either st, tu or us), 
and {och, • • • , /ch} are linear combinations of the (see Appendix B). 

For more general theories involving, e.g., fermions and/or derivative couplings, the nu- 
merator of Eq. (^ will also contain momenta in general. Work is in progress to extend 
the formalism to include cases, such as derivative coupling and fermions, where momentum 
dependence exists in the numerator. 

To illustrate our approach we will present here results for three choices of kernel: 

(a) Scalar-scalar ladder model with massive scalar exchange: The simple t-channel one-a- 
exchange kernel is given by 

Ap, q\ P) = 2 ■ (11) 

— [p — q) — 

The BS equation with this kernel together with the perturbative constituent particle prop- 
agator is often referred to as the "scalar-scalar ladder model" p|. Note that the kernel 
weight function is proportional to for this simple case. 

(b) Dressed ladder model: In this instance we dress the propagator of the exchanged a of 
case (a). The kernel then consists of the pole term as above, plus a piece which involves an 
integration (starting at a threshold of 4m^) over a mass parameter 7. 

(c) Generalised kernel: A sum of the one-a-exchange kernel Eq. (|TTp and a generalised ker- 
nel with fixed kernel parameter sets {7^*-*, ^''*''}. After the Wick rotation this kernel becomes 
complex due to the p ■ P and q ■ P terms, so that solving the BS amplitude as a function of 
Euclidean relative momentum would be very difficult in this case. 

The Wick rotated BS equation for cases (a) and (b) has been studied numerically 0. 
We use the previous numerical results for these kernels as a check of our new technique and 
numerical calculations. 
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B. PTIR for BS vertex 



As in Ref. lU, we will use the following form of the s-wave = 0) BS vertex: 

oo 1 

T{p, P)= [daf dz^^ , /"^"'"^i .,n , (12) 



-1 



m? + a -- [p^ + zp ■ P + ^] — ie 



with the following boundary condition for the weight function: 

hm = , (13) 

required to render the a-integral finite. A partial integration of Eq.(|T2|) with respect to a 
will serve to demonstrate that the positive integer is a dummy parameter, since it follows 
from such a step that weight functions associated with successive values of n are connected 
by the following relation: 

a 

p„+i(a,z) = n I da'pn{a',z) . (14) 





Note that the larger the dummy integer parameter n, the smoother the corresponding weight 
function. This is a particularly useful observation since we have a numerical solution of the 
BS equation in mind. 

Using the same arguments as outlined in it can be shown that the vertex PTIR for 
bound states with non-zero angular momentum i in an arbitrary frame of reference is 

°° 1 [(]( \ 

r[^A](p^ P) = y^p') [daf dz^^ 2! "'p^p2M, -in , (15) 

J [m^ + a — [p'' + zp ■ P + P^/4) — ie\ 

where P is an arbitrary timelike 4- vector with P^ = and p' = A~^{P)p. The Lorentz 
transformation A(P) connects P and the bound-state rest frame 4-vector P' = (M, 0), i.e, 
P = A(P)P'. The quantity yl^'ip') is the solid harmonic of order £, and may be written in 
the form \p'\^Y[^ {p'), where F/^ is the ordinary spherical harmonic of order i, p' = p'/\p'\, and 
where p' is the 3-vector relative momentum in the bound state rest frame. It is relatively 
straightforward to appreciate why Eq. (p!5D must be the correct form for a scalar bound 
state with angular momentum i. It follows from the self-reproducing property of the solid 
harmonics (see Eq. ( p6| ) in Appendix C) and from the fact that in the bound state rest 
frame p' is the only available three-vector. 

In the following sections we will study the BS equation Eq. (|l]) in an arbitrary frame 
in terms of this integral representation. Note that the dummy parameter n can always be 
taken sufficiently large such that the loop-momentum integral of the BS equation, Eq. (P, 
converges for any £ for which a bound state exists. 

Before proceeding, we should address the issue of so-called ghost states. A bound state 
whose BS amplitude ^{p,P) and equivalent vertex r{p,P) are anti-symmetric under the 
transformation P ■ p —P ■ p with fixed p"^ and P^ has a negative norm and is called a 
"ghost" 0. This symmetry corresponds to the one: z —>■ —z in PTIR form. We do not 
consider such states herein, as they are unphysical. 
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III. BS EQUATION FOR THE WEIGHT FUNCTION 



In this section we will reformulate the BS equation Eq. (|Tp as an integral equation in 
terms of the weight functions. This is the central result of this paper. We will very briefly 
describe the procedure, and state our main results; the details of the derivation may be 
found in Appendix C, as may definitions of the kernel and associated functions. 

We proceed by combining, using Feynman parametrisation, the integral representations 
for the scattering kernel and vertex with the bare propagators for the constituent particles. 
The procedure can easily be generalised to include dressed constituent propagators if desired, 
but we do not exercise this option here, for the sake of simplicity. After using the PTIR 
representations for the BS kernel and vertex {i.e., Eqs. (0) and (|15D ) in the right-hand 
side of the BS equation (Eq. (|I])), and after performing Feynman parametrisation for the 
right-hand side, the BS equation can be written as 



oo 1 



r[^/.](p,P) = y^^ f da I d-z _ \ 



oo 1 



X I da I dz^^^^^^ a" [A '°X^^\a, z; a, z)], (16) 



-1 



where F{a, z; p, P) = m? + a — {q^ + zp ■ P + jP"^) is a convenient shorthand notation. We 
have defined an "eigenvalue" A = g"^ / (47r)^, which we will use in our numerical work (see 
Appendix A). This has simply been factored out of the scattering kernel for convenience and 
for ease of comparison with other calculations in the ladder limit. The total kernel function 
'/C,{fl is defined in Appendix C (see Eqs. ( |C21| ) and ( |(J22| )) and its structure is discussed in 



tot. 



Appendix D with particualr attention paid to any potential singularities. 

Comparing Eq. (|TB|) with Eq. (^), and using the uniqueness theorem of PTIR 0, we 
finally obtain the following integral equation for p[f'(a;, z): 

= Ida /rfz-^M(a,^;a,.)^S^. (17) 

A a" J J ' a'^ ^ ^ 

-1 



This equation is the central result of this work. Note that in Eq.(|T^) we are solving for 
plfl/a""; this is for reasons of convenience for our numerical treatment of the BS equation. 

Summary: Since the weight functions PchiliC) for the scattering kernel are real functions 
by their construction, the total kernel function *°*/C|^1(q;, z] a, z) is real, so that Eq. (0) is a 
real integral equation in two variables a and z. Thus we have transformed the BS equation, 
which is a singular integral equation of complex distributions, into a real integral equation 
which is frame-independent. Once one solves for the BS vertex weight function, the BS 
vertex and the BS amplitude can be written down in an arbitrary frame. This is clearly 
advantageous for applications of the BS amplitude to relativistic problems. 
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IV. NUMERICAL RESULTS 



In this section we present numerical solutions for the BS vertex for bound states in scalar 
theories using Eq.(p!7D for three simple choices of scattering kernel: 

(a) pure ladder kernel with massive scalar exchange, 

(b) dressed ladder kernel with pole term as in (a), and 

(c) a generalised kernel combined with the pure ladder kernel of (a). 

The scattering kernel (a), i.e., the one-a-exchange kernel depending only on t = {p — qY, 
is given by Eg. ([TI|) . This corresponds to choosing for the kernel in Eq. (H) say: ptu = Pus = 
and in the st-channel 7 = ml, and a^j = c^t = 1, b^t = -2, 4t = e^t = fst = 0, [c.f., Eq. (0)], 
which amounts to choosing pst to be an appropriate product of 5-functions multiplied by 
g"^. In the pure ladder case it is convenient (and traditional) to factorise out the coupling 
constant and a factor of (47r)^, by defining the "eigenvalue" A = g'^/^AnY Thus it 

is usual to fix the bound state mass and then to solve for the coupling g^, which is what 
we have done here. 

In general (and for example, for the dressed ladder kernel) the kernel depends on higher 
powers of the coupling than g"^; in such instances we update the running value of A in an 
appropriate way during the iteration process (see Appendix A). 

Since the BSE is a homogeneous integral equation and we are only interested at present in 
extracting the coupling at fixed bound state mass, the choice of normalisation is unimportant 
provided that it is fixed in some reasonable way. 

The numerical solution of the vertex BSE is performed by choosing a suitable grid of 
a and z values, making an initial guess for the vertex weight function, and then iterating 
the integral equation (p!7| ) to covergence. One subtle point is that integrable square-root 
singularities may occur and must be appropriately handled numerically (see Appendices A 
and D for details). By optimising the choice of grid and increasing the number of grid points 
it was straightforward to increase the accuracy of the solutions to a relative error of 1 part in 
10^ and beyond. Further details of the numerical procedure used are given in Appendix A. 

A. Pure Ladder Kernel 

We have solved the vertex BS equation, Eq. ([T7|) , for a number of bound state masses 
between = (Goldstone-like bound state) and P^ = Am? (the stability threshold). 
Solutions were obtained for orbital excitations up to £ = 4 with no difficulties. We plot 
some examples of our solutions for £ = and £ = 1 in Fig. ^, and tabulate our results for 
the "eigenvalue" A = (7^/(47r)^ in Table |. A plot of the spectrum for £ = 0, i.e., of A vs. the 
fraction of binding rj = ^/P^/2m, is given in Fig. ^. All solutions presented were obtained 
using a a mass of ma = m/2. We have compared our eigenvalues to those obtained in the 
Wick-rotated treatment of Linden and Mitter , and have found agreement to better than 
0.03% for moderate choices of the {a, z)-grid. This is an improvement in accuracy of at least 
one order of magnitude over the results we have obtained previously for the BS amplitude. 
Furthermore, much greater accuracy is possible through an increase in the number of grid 
points used in the numerical integration should it be desired for whatever reason. 
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B. Dressed Ladder Kernel 



For the "dressed ladder" case, the scattering kernel is given by / = —iK = {ig)D„{ig), 
where is the renormalised cr-propagator at one-loop order. This is simply the sum of the 
pole term, Eq.(|lTD, and a continuum part, and is given by 



DAiP -<!)']= jds 



Pais) 



s — [p — qY — ie' 



(18) 



where Po-(s)=ppoio(s) + Pcontls). Ppoie(s) is simply 5{s — rri^), and it can be shown that, to 
one-loop order, pcont{s) is given by the following expression: 



Pcont = X9{s- Am? 



' s — Am? 1 



Ms) 



(19) 



where A(s) is the following function: 



A(.) 




Am^ — mi 



arctan 



mt 



mt 



m 



y^m^(4m^ — 



^ 4m2 - 
arctan 



mt 



mi 



^ 4m^ 



mi 



2„2 



Is — Am^ , / + ys — 477? 
+ 2W In 



s — Am^ 



(20) 



Note that the use of Eq. (^) introduces an extra integration (over the mass parameter 
s). We performed this numerically using gaussian quadrature; 10 to 15 quadrature points 
in s provide solutions of satisfactory accuracy. 

We have solved Eq. ([T7|) for the dressed ladder kernel above, with the pole being situated 
at ttIq. = 777, for various values of the bound state mass squared, P^. As for the pure ladder 
case above, solutions have been obtained up to £ = 4. For example, our s-wave {C. = 0) 
eigenvalue for = 3.24m^ and for the exchange particle pole at m^ = m is X = 1.516 
for a grid of 80 x 41 and is A = 1.518 for a grid of 150 x 91. The corresponding Linden 
and Mitter value is A^; = 1.518. Even with a relatively coarse grid high accuracies result. 
Similarly, we have found for other values of P^ that an accuracy of 0.3% or better is routinely 
attained, even with the use of the coarse 80 x 41 grid. As the above results demonstrate 
higher accuracy is easily obtained at the cost of more CPU time. 

We plot the re-scaled weight function, p§^/a'^, for the case P^ = 0.04777^ in Fig. | for 
purposes of comparison with the corresponding ladder case. We have chosen this value of P^ 
since a smaller bound-state mass corresponds to tighter binding and hence larger coupling, 
which should enhance the effect of this one-loop self-energy insertion. We see that the shape 
of the weight function is not qualitatively very different from the ladder case. 
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C. Generalised Kernel 



This particular example of a "generalised" kernel is an instance of a scattering kernel 
for which Euclidean space solution is not possible. We have solved the BS equation for this 
case, in particular for a sum of the pure ladder kernel as described above and two randomly 
chosen non-ladder terms, each with a weight of 0.25A, and with the following fixed parameter 
sets (in the st channel): 

sty ^st) dst, ^st, fst} — 

(1) {2.25m2, 0.47261150181, -0.29743163287, 0.58277042955, 
0.28282145969, -0.23965580016, 0.32196629047} 

(2) {2.25m2, 0.47261150181, 0.29743163287, 0.58277042955, 
0.28282145969, -0.23965580016, -0.32196629047}. (21) 

The parameters {a^t, . . . , fst} were obtained from a set of values . . . , ^e} produced by 
a random number generator (see Appendix B). This was done to emphasise that our tech- 
nique produces well-behaved solutions for an arbitrary kernel. We have had little difficulty 
obtaining noise-free solutions for this kernel for orbital excitations up to £ = 4. 

This kernel yields an s-wave eigenvalue of A = g'^/{47r)'^ = 1.3569 for = 1.44m^, 
c.f. the ladder value of A = 1.9402 for the same bound state mass. We therefore find, as 
before |jl|, that the additions to the pure ladder kernel have enhanced the binding, i.e. they 
are attractive. This is of course to be expected in a scalar theory and was also observed 



in Refs. [|TO|,|TT|. Not only is the eigenvalue lower, but additional structure is present in 
the vertex weight function (see Fig. In contrast to the solutions obtained for the BS 
amplitude in the previous work [|l|, we find that there is no numerical noise in the vertex 
weight function. If one compares the general kernel example solution in Fig. ^ with the 
ladder solution for the same bound-state mass (see Fig. 0), it is readily apparent that there 
is some additional structure superimposed upon the weight function, due to the addition of 
the generalised kernel terms. 



V. SUMMARY AND CONCLUSIONS 

We have derived a real integral equation for the weight function of the scalar-scalar 
Bethe-Salpeter (BS) vertex from the BS equation for scalar theories without derivative 
coupling. This was achieved using the perturbation theory integral representation (PTIR), 
which is an extension of the spectral representation for two-point Green's functions, for both 
the scattering kernel [Eq. (^j and the BS vertex itself [Eq. ([T5|)]. The uniqueness theorem 
of the PTIR and the appropriate application of Feynman parametrisation then led to the 
central result of the paper given in Eq. (|1^). 

We have demonstrated that Eq.([T7|) is numerically tractable for several simple kernels, 
including a randomly chosen case where it is not possible to write the kernel as a sum of 
ordinary Feynman diagrams. Our results for both the pure and dressed ladder kernels are 
in excellent agreement with the results obtained previously in the Wick-rotated approach. 
The agreement for the pure ladder kernel is even better than in Ref. |]l[, vindicating our 
decision to solve here the vertex equation rather than the amplitude equation. We obtained 
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an accuracy in all of our results of approximately 1 in 10^ with modest {a, z)-grid choices on 
a workstation. This can be improved by using finer grids and larger computers as desired. 

Further applications of our formalism are currently being investigated, particularly the 
crossed ladder and separable kernels. These will not only provide yet another test of our 
implementation of the method, as Euclidean space results are also available for these cases, 
but also will provide us with an opportunity to solve a problem featuring more realistic 
scattering kernels. 

It is also important to consider how the PTIR can be extended to include fermions 
and derivative coupling, so that we have a covariant framework within which to study, 
for example, mesons in QCD using a coupled Bethe-Salpeter — Dyson- Schwinger equation 
approach. This would require us to incorporate confinement into the PTIR, which at this 
stage remains another important and interesting challenge. 

[Those interested in applications of this technique may request a copy of the 
computer code from AGW at the given e-mail address.] 
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APPENDIX A: ALGORITHM 

Here we detail the algorithm used in our numerical studies of the integral equation: 

1 ^ = / %i) «kL'>(^.. a, A). (Al) 

-1 

where we have explicitly shown the coupling dependence of the kernel function. We have 
defined here an "eigenvalue" A = (7^/(47r)^. Our rationale for this is as follows: for a 
given scattering kernel the integral equation ( |A1| ) may be solved for the bound state mass 



= M^. However, the dependence of the kernel function *°*/C|f](Q;, a, 2;) on the bound 
state mass is highly non-linear and complicated. It is therefore convenient and traditional 
to instead solve the equation for the coupling (7^, which appears in the weight function 
Pch{l,0 ^he scattering kernel, with a fixed bound state mass P^. We first fix the bound 
state mass P^ and regard the integral equation ( [Al|) as an "eigenvalue" problem. The 
"eigenvalue" is then introduced by factorising the coupling constant from the scattering 
kernel weight function Pch(7,0- this convention, the kernel function ^"%]^^(a,z;a, z; X) 
becomes a power series in A starting from 0{1) for a perturbative scattering kernel. 

Strictly speaking, the integral equation ([Al| ) is not an "eigenvalue" equation, since the 
kernel function itself contains A in the general {i.e., non-ladder) case. We thus solve the 
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equation by iteration rather than applying methods for eigensystems. With an appropriate 
initial guess for the weight function p|f'(a, z) for the BS vertex and the coupling constant A we 
generate the new weight function by evaluating the RHS of the integral equation (|Al|). The 



"eigenvalue" A associated with the weight function is extracted by imposing an appropriate 
normalization condition which we will discuss later. This generated weight function and its 
"eigenvalue" are used as inputs and we obtain updated values, which ought to be closer to 
the solution than the input values, by evaluating the integral. We repeat this cycle until 
both the "eigenvalue" and the weight function converge. 

The normalization condition for the BS vertex or equivalent BS amplitude in momentum 
space is well known and involves the derivative of the scattering kernel with respect to the 
bound state 4-momentum P (See Eq.(^ in Sec.|TD. When the scattering kernel depends 
on the total momentum P, the normalization condition is the integral over two relative 
momenta; one for the BS vertex and the other for the conjugate one. The corresponding 
normalization condition in PTIR form is written as the 4-dimensional integral over spectral 
parameters a and z. Imposing a condition that involves such a multi-dimensional integral in 
the iteration cycle makes the calculation less accurate and time consuming. We shall rather 
use a suitable normalization condition for the weight function p,[f^(a;, z) during the iteration. 
The physical normalization condition may be imposed by appropriately rescaling the 
obtained solution. Of course, the value of A = g^/^An)^ is unaffected by the choice of 
normalisation of the vertex weight function. 

Since we are considering bound states whose constituents are of equal mass m, we ex- 
pect that a physically reasonable scattering kernel I{p, q; P) will give a kernel function 
/C|f](a, z] a, z] A) symmetric under the transformation z — > —z. The weight function pif^(a, z) 
is then either symmetric or anti-symmetric m. z ^ —z. For a symmetric solution the fol- 
lowing normalization is convenient: 

da I dz ^£1^ = 1, (A2) 



provided that the integral does not identically vanish. From Eq. ( |C21D , the kernel function 



is given by the difference of two terms, viz /C|fl(a, 2;; 0, 0; A) — /C|f'(a, z; a, z; A). With the 
normalization condition ( |A2| ) the integral equation ( |A1| ) can be written as an inhomogeneous 



one: 



00 1 



/CM(a, z- 0, 0-X)- fda f dz /cM(a, z; a, z; A) 



1 

A J J a 



(A3) 



On the other hand, the integral of the weight function over a and z for an anti-symmetric 
solution vanishes identically. Then the integral equation ( |A1| ) is written as 



A a" J J a 

-1 



As discussed in the main text, the bound state whose BS amplitude $(p, P) and equivalent 
vertex T{p, P) are anti-symmetric under the transformation P ■ p —>■ —P ■ p with fixed p^ 
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and has a negative norm and is called a "ghost" This symmetry corresponds to the 
one: z —z in PTIR form. Thus a bound state whose weight function is anti-symmetric in 
z-reflection and which satisfies the homogeneous integral equation ( [A4| ) is a "ghost" state. 
We hereafter concentrate on the "normal" solutions, namely 2;-symmetric ones, and on the 
integral equation ([A3|) . 

It can be shown that the inhomogeneous term vanish unless a > ath, where ath is the 
threshold of the weight function depending on the value of z for a given scattering kernel. 
For a one-(T exchange kernel with the mass fi, the threshold can be written as 



m 



1/2 



m'-{l-z') — 



(A5) 



This threshold determines the support of the weight function p or equivalently (f for the 
normal solution. Although we cannot write the threshold in a simple form like ([A5|) for 
general scattering kernels, we can extract it numerically by analyzing the inhomogeneous 
term. On the other hand, the kernel function /C|fl(Q;, z; a, z; A) has the support property for 
a given a, z and z that it vanishes unless a is less than some value Omax- For the case of a 
one-o" exchange kernel it is given by 



^ma,x\CX, Zj 



a + ( - (1 - 



P 



2\ 1/2 



/i 




As in the case of amin the analytic form of the upper limit Omax is unknown, so we extract the 
corresponding upper limit numerically for general scattering kernels. Writing these limits of 
the integral explicitly the integral equation we use is then 



1 P[fl(«,^ 



A a' 



K,l\a, z- \) - dz da ^" \' ' fCl\a,z-a,z-\). (A7) 

J J at ' " " ' ' 

-1 



a'- 



We evaluate the integral over a and z in the RHS of Eq. (^3) follows. Recall that the 
kernel function K\^{a, z\ a, z; A) is given by the following integral: 



K^nia, z; a,z;X)=J2 J j d-f 

n 



— Pch(7,6 
9 



Kl^\a,z;a,z;-f,^) 



(A8) 



We thus start by replacing the integrations over Feynman parameters and the spectral 
variable 7 by summations over discretised variables. Secondly, we map the semidefinite 
range of a G [0, 00) to the finite one y G [0, 1]: 



a = ao + C 



y 



1-y 



(A9) 



where ao and C are some constants which should be chosen such that the weight function 
is largest around the mapped variable y ~ 1/2. We then discretise both y (equivalently a) 
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and z and prepare the initial weight function on this grid. For each cycle of the iteration 
we perform the integral as follows. We first evaluate the a integral for a given point in 
the a and z plane and on the z grid. For each value of the discretised ^ and 7 we extract 
the support of the the kernel function -ft'lf' ( A). We then divide the integral range 

[ath(2;), Q;max(Q;, ^, 2)] into subranges according to the support of the kernel function with 
discretised ^ and 7. As discussed in Appendix D the kernel function ii"]fl(a, z] a, z] A) may 
diverge as an integrable square root singularity at the boundary of the support. While 
this is always the case for the one a exchange kernel, the kernel may take a finite value 
in general. We thus choose an appropriate integration method to perform the integration 
over a for each subrange. The weight function at arbitrary a is evaluated by interpolating 
the values of p on the grid. We perform the a integral in this way for each grid point 
of z and the integral over z is performed by interpolating these values. With this careful 
treatment of integrable square root singularities we need not introduce any regularization 
or cutoff parameters. Furthermore, this method allows us to choose the a and z grid for the 
newly generated weight function independent of the a and z grid. We optimise the "new" 
grid by analyzing the shape of the "old" weight function used in the RHS of Eq.(^). The 
eigenvalue is evaluated using the normalization condition ([A2|) . 



APPENDIX B: PTIR FOR SCATTERING KERNEL 

In this appendix we list the dimensionless coefficients {och, &ch, Cch, • • • , /ch} in Eq. (p!OD 
for different channels {ch} = {st}, {tu}, {us} in terms of the Feynman parameters defined 
in Eq. m. 





St 


tu 


us 




6 + 6 + 6 


6 + 6 + 6 + 6 


6 + 6 + 6 


hch 




2(6-6) 


26 


Cch 


6 + 6 + 6 


6 + 6 + 6 + 6 


6 + 6 + 6 


C^ch 


1(6 + 6 + 6 + 6) +6 


1(6 + 6 + 6 + 6) 


i(6 + 6 + 6 + 6) + 6 


Cch 


6-6 


6-6 


6-6 


fch 


6-6 


6-6 


6-6 



Beginning with the above definitions for the scattering kernel parameters in terms of the 
Feynman parameters ^, and noting that I]^=i 6 = 1; it is possible to prove the following 
relations between the kernel parameters, for all three channels: 



ac > 

4 

2 



< c 



af 



eb 
~2 



< ac 



4 ' 



(Bl) 
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APPENDIX C: KERNEL FUNCTION 



In this appendix we detail our derivation of the real integral equation for the BS vertex. 
We begin with the following PTIR form of the bound-state vertex [0: 



oo 1 



r[^'^^](g,P) = y',^{A-\P)q) jdajdz 



-1 



[F{a,z]q,P) - ie] 



(CI) 



In Eq. (pl|), yl^"{Ar^{P)q) is the solid harmonic for a bound state with angular momentum 
quantum numbers £ and ^z■, ^) is the PTIR weight function for the bound state vertex 

function, and n is a dummy parameter. The function F is given by 



z]q,P) = a + 



^^'^m'-iq+ \Pf) + - (-g + \pf) 



2 2- / / ■ 2 

a + m^ - {q^ + zq- P + ip^). 



(C2) 



We proceed by substituting this form of the vertex into the vertex BSE, Eq. (P, and 
combining the various factors on the right hand side of the resultant equation (bare propaga- 
tors, scattering kernel and vertex PTIR) using Feynman parametrisation. We first combine 
the bare propagators for the scalar constituents with the denominator of the vertex PTIR: 



D{q+\P)D{~q+\P) 



1 T{n + 2) 
2r(n)r(2) 



[F(a,z; g,P)-ze]« 
1 1 

j dri j dtr-\i-t) 



-1 



[F{ta, tz+{l- t)r]; q, P) - ie]"+2 ' 



(C3) 



We now combine the factor 1/[P(. • .)]"''~^ from the integrand with the denominator of 
the PTIR for the scattering kernel (see Eq. (0)). After some algebra this yields 



1 



7 - (ag2 + bp-q + cp'^ + dP'^ + eq ■ P + fp ■ P) - ie [F{ta, tz + {1 - t)r]; q, P) - ieY+'^ 
r(n + 3) ' 



r(i)r(n + 2) 



dx 



X 



n+l 



1 -xY+^ (|/-Fa)"+3 



1 



[{[c{y + a) - {h^lA)]/{y + af]F' -{q+ [bp + (e + [tz + (1 - t)7]]z)P]/2{a + y)y - te]^+^ ' 
where 

F' = F[A{ta, tz + {l- t)7]; y), Z{tz + (1 - t)r/; y); p, P], 
1 



A{a,z] y) 



c{y + a)-'^ 



a + m" -{l-z')^\{y + af 



P2 



7 — (a + c)m — aa + {a + c — Ad + 2z{e — az))— > {y + a) 



+ — Im^ - \1 



-2(az — e) 



2\ p2 
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c{y + a) - 



4 



Z/ = T^- (C5) 

1 — X 

Note that the parameters {a, ...,/} do not exphcitly have the subscript "ch" attached 
to them in this instance (cf. Appendix B and Eq. (^), for the sake of brevity. We will also, 
for the time being, omit for brevity the sum over channels st, tu, us in the kernel function. 

Having combined all factors on the right hand side, we are now in a position to perform 
the integral over the loop momentum q. To do so, we must utilise the following property of 
the solid harmonics, 

/ d'qHf)yHq + p) = yHv) J d'gnf), (ce) 

where F is a sufficiently rapidly decreasing function which gives a finite integral . Also note 
that A-^(P) boosts the 4-vector P to rest, i.e. A-^{P)P = {y/P^,0), and that the solid 
harmonics are functions purely of the three- vector part of their argument, i.e. y^^ip) = 
\p\^Yg''{p) (here p = p/lpl). Bearing this in mind, we obtain for the loop integral 



/ 



(2n)H 



n+3 



( hV 1 r(n + i) yl^{k-\p)p) 

''\ 2) (a + y)^ (47r)2r(n + 3) (M2 _ ' 



(C7) 



where is simply that part of the denominator of Eq. (|C^) that does not depend at all 
on the 4-momentum q. 

Ignoring integrations over weight functions for the moment, we have, after performing 
the loop momentum integral, the following result: 



I 72^^^^' ^' ^^^^^ + \P)D{~q + \P) 



(27r)42 ' ' 2 ^ V ^ ' 2 '[F{a,z-q,P)-ie] 

1 1 / hYvin + l) 



(47r)2 2 V 2/ r(n 
X I dr] I dte-\l-t) / dx 



11 1 , 1 / N 1 

1 yu+l ^ 



'10 



1 _ a;)2 + ay [c{y + a) 



_ ^In+l 



yH^-\p)p) .C8) 

[F{A{ta, tz + {l- t)v; y), Z{tz + (1 - t)r^- y)-p, P) - lef^' ' 

In order to obtain a real integral equation involving only weight functions, it is necessary 
to recast the last factor in Eq. ( |C8|) in a form similar to that found in the vertex PTIR, 
Eq. (pip. To proceed we therefore insert the trivial integral 

1 

J dz S{z- Z{tz + {l-t)r];y)) = 1 (C9) 
-1 
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into the right-hand side of (|C8|) , and ehminate the integration over t by rewriting the delta 
function in terms of t. We are permitted to do this because the function Z is bounded 
between -1 and 1. That this is true is easily seen by observing that Z is monotonic in the 
variable y, with y G [0, oo), and then by taking the limits y ^ and y oo. The former 
limit gives 



Z{z;y^O) 



(CIO) 



and so from the third inequality in Eq. ( [Bl| ) we have that \Z{z;y — ^ 0)| < 1. The second 
limit gives 



2^iz;y 



oo 



(Cll) 



which allows us to use the second inequality in Eq. ( pT]) to conclude that \Z{z; y oo)\ < 1, 
given that ( |C11|) is monotonic in z. Since \Z\ < 1 in these two limits, Z must be bounded 
between —1 and 1 for all y G [0, oo). 
The insertion of this integral gives us 

1 1 
dz 



1 



(47r)2 2 



T(n + 1] 



T(n) 



dx 



,n+l 



(y + a 



(1 - a;)2 {y + aY h^+^ 



G 



n-1 



X {e{z-G)e{G + l) / dr] 



h z-G fG-ri^ 
y{z-riY \ z-ri ) 

I y{z-riy \ z-ri ) 
where we have introduced the following: 

G-T] 



[FiAitoa,G),z;p,P)-te]--^^ 
[F{A{toa,G),z;p,P)-te]-+^ 



tn = 



Z — 7] 



G^G{z-y) 



1 



-2{cz-f)^ , , b 

^ (l/ + a)-(-^^ + e) 



h = h{y) = c{y + a) - —. 
We next make a change of variable rj —>■ a, such that 



a 



yiy + a). , 9{z]y) 

-toa + 



h ^ h ' 

g{z- y) ^ A{z){y + af + B\z){y + a) + G'{z) 



^{z) = m 

B\z] 
G\z\ 



1 



'-2{c-z-jy 



p2 

T 



7 - 
62 



(a + c)vr? + (a + c - 4(i - 2z[cz - f))-^ 



-im 



p2 

T 



(C12) 



(C13) 



(C14) 
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Note that for brevity we sometimes write g{z; y) as g below. This should not be confused 
with the coupling strength g, since the meaning should be clear from the context. 
With this transformation the factor in braces in Eq. ( |C12D becomes 



Rj^a+g/h 



e{z-G)9{G + l) J 



da 



1 



g/h 

R^a+g/h 



-d{G - z)e{l - G) J 



da 



g/h 



a'^ y {y{y + a)y 
1 1 

a" y iyiy + a)) 



iha-gr-' 
-{ha-gr-' 



X 



[F(a,^;p,P) -ie]"+i' 
with the functions R± being defined by 



yiy + a)Giz;y)±l 



h{y) z±l 

We may use the following to shift the limits of integration of a to (— oo, oo) 

b OD oo oo 

J da = J da — J da = J da[9{a — a) — 6{a — b)]. 



(C15) 



(C16) 



(C17) 



The expression ( P15|) then becomes 



Y.O{±iG-zmiTG) [da^- 

^ —no 



- 9 



n-l 



X 



eia- 



9 a 



9_ 

h 



R^a 



[F{a,z;p,P) - ie]"+i' 



(C18) 



We complete our derivation of the integral equation by integrating by parts with respect 
to a in order to reduce the power of 1/F{. . .) from n + 1 to n, noting that the boundary 
term resultant from such an integration vanishes due to the presence of the step functions. 
We therefore have, finally, 



oo 1 

r[^'^^i(p,p) = 3;f (A-i(p)p) Jdajdz 



1 



[F{a,z;p,P) -ie] 



X da dz 



a 



a' 



A *°*/CM(a,z;a,z) 



(C19) 



We may use the uniqueness theorem of PTIR to obtain the equation which we will 
solve numerically: 



-1 



a' 



(C20) 
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where the full analytical expression for the kernel function *°*/C can be written as 

*°*/C = /CM(a, z] 0, 0) - K,l\a, z] a, z), (C21) 
where K, is the function: 



00 

^; a, ^) = H / di j d-f 



* Q 



|6| \ 2j dal ' {y'y+^ [ h'{y') J 
xJ2e{±{G{-z;y')-z))e{lTG'{-z;y')) 

( g'iz,y) y'{y'-a)lTG'{z;y') \ 

Note that we have made a shift of variable from y to y'. The quantities ^f', h' and G' are 
the same as their unprimed counterparts, except that the y dependence of these functions 
has been transformed according to y ^ y' — y + a. Note also that we indicate explicitly the 
dependence of Klf^ on the scattering kernel parameters {7, C}- For the remainder we will 
omit these additional labels for brevity. 

In order to implement this kernel numerically, we must perform the derivative with 
respect to a and simplify the resultant expression, as well as transforming those integration 
variables with semi-infinite or infinite ranges to variables which have a finite range. We 
begin by performing the a derivative, which splits the kernel into two pieces, one of which 
contains a delta function. After this differentiation, we have 



"2) 



The piece containing the delta function may be integrated over y' in a relatively straight- 
forward manner, simply by rewriting the delta function in terms of y'. The argument of the 
delta function is quadratic in y': 

argument = 7o-{Azf{z; a, z)y''^ + Bzf{a, z; a, z)y' + C(q;, z)). (C24) 

The delta function is therefore 

5 ( ^{A^y" + B^y' + C)\=jz I A y' - y'MD^{a, z; a, z)), (C25) 

V ~ T / 1=1 y D^{a,z;a,z) 
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where = — ^A^C, and the y[ are the roots of the quadratic, i. e. 



2A^ 



y2 = 



(C26) 



The case n 2 is of particular interest to us, and so we will restrict ourselves to this 
case from now on. Dropping the prime on y, the n — 2 kernel function may be written as 



\o\ 



1 1 



X 



cy-'i 



(^A^(z; a, z)y^ + B^:(a, z; a, z)y + C(q;, . 



^ aO {D^{a, z] a, z)) ^ 



It 



-2(cg-/) 



aT -f^ + e 1 



^jD^{a, z\ a, z) i=i\ 1 T 



l^Z 



Vi 



dy6{y - yi) 



X e 
X e 



1 



± 



-2(cz-/) 



- z\y±{az - I "2^ + ^ 



1 T — ^ — — \ y-a± \ --Z + e 



(C27) 



where 



Dzf{a^ z\ a, z) = B-^{a, z; a, z) — AA^:{z] a, z)C(a, z) 



A^{z] a, z) — — 1 



-2{c-z-fy 



P2 1^ 



-2(cg-/) 



l^Z 



-a 



Bzf{a, z;a,z) = 7 — ca — (a + c)m + (a + c — 4d — 2z{cz — f))— 

C{a,z)^j{a + m'-{l-f)^). 



p2 (-|^ + e) 



l^z 



a 



(C28) 



For the purposes of numerical solution we now make successive transformations of the 
integration variable y, first to y = 1/y, and then from y to Y = (6^/4)^. The first trans- 
formation serves to render the range of integration finite, while the second ensures that we 
do not encounter any difficulties in the kernel function in the limit 6^0, which can occur 
for example in the separable kernel case. The kernel function after these transformations 
becomes 



+ 



Dzp{a, z;a, z) i=i 



YfJdY5{Y-Yi) 
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X e {±{az - ho)Y ± (^0 - j^)) ({-a ± ho)Y + j T ^o) , (C29) 

where ho{z) = {—b/2)z + e, and ^o(^) = i~b/2){cz — f). The Yi are the roots of the quadratic 
C^Y^ + B^Y + A^, and 

C^{a, z; z) = {lT z) (^a + - {1 - ^^)-^^ 

~ _ _ _ _ _ p2 

-B=p(a, z; a, 2;) = (1 =]= z){'j — ca — {a + c)m? + {a + c — Ad — 2z{cz — /))— ) 

- =F (-^^ + e))a 

A^{z- a, z) = {lT z) [^-{m' - — ) + {cz - ff—^ + ^ T^'" " " 

D^ = Bl- AA^C^ (C30) 
This is the expression which we implement numerically. Note that the support of the kernel 



is entirely determined by the step functions in Eq. (|C29| ) . In general it is not possible to 
extract the support analytically, and so in most cases this step must be done numerically. 



APPENDIX D: KERNEL SINGULARITIES 

In this section we discuss the structure of the kernel function K]l\a, z\ a, 2;) for arbitrary 
£ with a fixed kernel parameter set (7,0; i-^-, for constant {7, ach, &ch, Cch, • • • , /ch}- We 
will in this section omit for brevity the subscript ch. Since the case n = 2 is of particular 
interest to us for numerical treatment, we discuss possible singularities of the kernel function 
K^n=2^'^-: z\ a, z), whose expression and derivation are given in Appendix C. General n cases 
can be also considered in a similar manner. 

As shown in Appendix C the kernel function K^}^2{'^^^'i'^jz) given by the Feynman 
parameter integral consists of two terms, one containing only step functions and another 
containing a delta function. It is convenient to make the Feynman parameter y finite to 
discuss the singularities of the kernel, and so we will discuss the structure of the kernel 
function based on the expression ( |C29|) . 

The step function term is given by the integral 

dYY^-\ (Dl) 

Y . 

^ mill 

where the upper Fmax and lower limits Y^^^ of the integral are determined by relatively 
complicated step functions depending on the variables a, z, a, and z as well as the kernel 
parameters {7, a,. . .,/}. It is easy to show that the upper limit Fmax is finite as long as 
the parameter a does not vanish. Since the scattering kernel with identically vanishing a is 
nothing but the constant scattering kernel in the relative momentum p, we do not consider 
this case. Thus the Feynman parameter integral may diverge logarithmically, and this only 



if Ymin vanishes for the i = case. However, as is clear from the expression (|C29| ), the point 
y = is always excluded by the step functions, so this integral never diverges. 
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The delta function term can be written as a sum of fractions with square root factors 
in their denominator together with finite numerators. The square root factor comes from 
the Jacobian to change the variable of the delta function from the spectral variable a to the 
Feynman parameter Y . Note that this situation is quite general and occurs for any angular 
momentum i and dummy parameter n. From the argument of the square root, this term 
becomes singular if a satisfies the following quadratic equation: 

-4C(a,.lM(.1 + + ^0, (D2) 

a=F/io(^) V a=F/io(^) / j 

where h^iz) = {—b/2)z + e, and go{z) = {—b/2){cz — /). The functions A{z), B{a,z) and 
C{a,z) are 

B(a, = 7 — ca — (a + 0)711^ + (a + c — Ad — 2z{cz — /)) — , 

p2 

C{a,z) = a + m^ - {1-z^)—. (D3) 

Thus the kernel function diverges as a square root if Eq. ( P^ ) possesses a simple root. On 
the other hand, the kernel function diverges linearly if Eq. (|D2| ) admits a double root. Since 
C{a,z) > for any bound state, a double root occurs only if the terms in the second set 
of parentheses cancel. In this case, however, the residue of this pole (linear singularity) 
vanishes, so that the delta function term stays finite as a whole. 

To summarise: the kernel function ii']fl2(Q;, ^; a, 2;) for a fixed kernel parameter set 
{7, ttch, bch, ■ ■ ■ , fch} contains only integrable square root singularities at the boundary of 
its support, which if appropriately treated numerically present no difficulties. 
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TABLES 



rj 


Xe 


A 


7] 


Xe 


A 





2.5658 


2.5662 


0.8 


1.4055 


1.4056 


0.2 


2.4984 


2.4988 


0.9 


1.0349 


1.0350 


0.4 


2.2933 


2.2937 


0.99 


0.5167 


0.5168 


0.6 


1.9398 


1.9402 


0.999 


0.3852 


0.3853 



TABLE L Comparison of the coupling strengths, A = g^/(47r)^, obtained for the ladder ap- 
proximation kernel from the Euclidean (i.e., Wick-rotated) s-wave solution (Xe) and those ob- 
tained here directly in Minkowski space (A), using a moderate grid choice for a and z. The 
Wick-rotated values are from Linden and Mitter . The parameter ry is the "fraction of binding" , 
rj = V P'^ /2m = M/2m, where = and M is the mass of the bound state (i.e., < ?7 < 1). 
The results shown here had = m/2 for the exchange particle (a) mass. 
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FIGURES 



iD 




iD 

iD 




iD 



FIG. 1. Diagrammatic representation of the Bethe-Salpeter equation for (a) the BS vertex (r) 
and (b) the BS amphtude {^). The fully dressed constituent particle propagator is denoted by D 
and K is the scattering kernel for the constituents. 
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(a) (b) 

FIG. 2. Solutions for the bound-state vertex weight function for the s-wave {i = 0) and p-wave 
= 1) cases are given in (a) and (b) respectively for a variety of values for the fraction of binding 
rj = ■\/P^/2m = M/2m. The exchange particle (a) mass is rriu = m/2 in these solutions. It is 
convenient to plot the rescaled weight function p{a,z) = P2\a, z)/a'^ in these figures. 
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Eigenvalue Spectrum 
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FIG. 3. The s-wave (i.e., ^ = 0) bound-state spectrum for the ladder kernel, where A = 5^/(47r)^ 
and r] = V P'^ /2m = M/2m. The exchanged particle (cr) mass is half the constituent mass, 
TTio- = m/2, in these results. These results correspond to those in Table |. 
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